Nonlinear spin Hall effect in GaAs (110) quantum wells 
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We consider stationary spin current in a (llO)-oriented GaAs-based symmetric quantum well 
due to a nonlinear response to an external periodic electric field. The model assumed includes the 
Dresselhaus spin-orbit interaction and the random Rashba spin-orbit coupling. The Dresselhaus 
term is uniform in the quantum well plane and gives rise to spin splitting of the electron band. The 
external electric field of frequency u) - in the presence of random Rashba coupling - leads to virtual 
spin-flip transitions between spin subbands, generating stationary pure spin current proportional to 
the square of the field amplitude. 

PACS numbers: 72.25.Rb, 72.25.Hg 



I. INTRODUCTION 

Spin-orbit (SO) interactions in semiconductors reveal a 
variety of fundamental spin related phenomenal Forma- 
tion of stable spin helices with nontrivial temporal and 
spatial dynamics^— spin optics^ and spin-dependent 
sound radiation^ represent only some of these effects. 
Spin currents* 7 - solely attributed to the spin-orbit cou- 
pling, provide a possibility of inducing and controlling 
spin motion by electrical and optical fields, and there- 
fore became one of the key elements of the modern spin- 
tronics oriented at new device applications of semicon- 
ductor based structures. Thorough investigations of re- 
alistic systems enlarge the variety of both fundamental 
phenomena and possible applications. As an example we 
mention that disorder - always present in real systems - 
plays a crucial role in the spin-Hall effect^, as the spin 
Hall conductivity can be totally suppressed by any finite 
concentration of impurities. 

Recently, symmetric GaAs (110) quantum wells (QWs) 
became the subject of extensive experimental and the- 
oretical investigations. This is related to the expec- 
tation of the longest spin relaxation times in these 
structures^— which in turn can lead to interesting spin 
dynamics i 17 ' 18 A stationary pure spin current accompa- 
nying an electric current in (110) QWs was observed as 
reported in Ref. [l9[. SO interaction in these systems, 
described by the Dresselhaus term in the corresponding 
Hamiltonian^Sr— conserves the electron spin along the 
axis normal to the QW plane for any electron momen- 
tum k. As a result, random motion of an electron does 
not lead to a random direction of the spin-orbit field 
and therefore does not lead to spin relaxation. In re- 
ality, however, this spin component relaxes very slowly, 
and its analysis provides a test for the rapidly developing 
low-frequency spin-noise spectroscopy^ suitable for the 



measurements of the long spin evolutions. 

In the case of perfect z — > — z symmetry (the axis z is 
perpendicular to the QW plane), the Rashba SO inter- 
action is zero. In real structures, however, the Rashba 
coupling still exists in the form of a spatially fluctuat- 
ing SO field (though being zero on average) 2^— This 
interaction induces spin-flip processes leading to the spin 
relaxation^ and can be also responsible for generation 
of a nonequilibrium spin density due to the absorption 
of an external electromagnetic fields Recently, it was 
proposed that this random SO coupling can result in the 
spin orientation by an external current^ 7 - and also can 
play a role in the formation of the stripe structure of 
spin current distribution^ 

In this paper we propose a new possibility of exciting 
a steady pure spin current by a periodic external filed, 
extending thus the abilities of spin manipulation in real 
situations. In contrast to the conventional spin Hall ef- 
fect, which is linear in the external electric field, the pro- 
posed spin current is quadratic in the external periodic 
field. The effect is a result of the interplay of constant 
Dresselhaus and spatially random Rashba terms, and is 
not related to the spin currents produced by gate ma- 
nipulation of the Rashba coupling^ or adiabatic pump- 
ing in graphene^S Exact mechanism of the effect does 
not necessarily involve real spin-flip transitions of elec- 
trons between the spin-split subbands in (llO)-oriented 
GaAs QW, but relies on virtual spin-flip processes which 
rcnormalize the wave functions of electrons in a nonequi- 
librium state. This makes such a nonlinear current a 
physically new phenomenon, which appears if one ac- 
counts for more realistic effects than those described by 
the conventional Rashba and Dresselhaus models. 

In Section 2 we describe the model and Hamiltonian 
of the system. Spin current is calculated in Section 3. 
Summary and final conclusions are presented in Section 
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II. MODEL 



aak x . The energy bands of spin up and spin down elec- 
trons are thus shifted in opposite directions along the k x 
axis. The corresponding Green function is then diagonal 
in the spin subspace, 



Hamiltonian of a two-dimensional electron gas with the 
constant Dresselhaus term Hjj and spatially fluctuating 
Rashba spin-orbit interaction Hr, subjected to external 
electromagnetic field described by the vector potential 
A(r, t), takes the following form (we use units with fi = 1) 



H = H + H L 



where the first two terms are 
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where a = + for spin up (t) electrons and a = — for spin 
down (J,) electrons, whereas <5k CT is the momentum and 
spin dependent relaxation rate. 
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The Dresselhaus constant a = r yw 2 /2w 2 , where 7 is the 
corresponding bulk Dresselhaus coupling parameter, is 
inversely proportional to the square of the QW width 
w. The other components of the Dresselhaus interac- 
tion vanish due to the specific symmetry of the (110) 
orientation . 16 ! 21 

The last term in Eq. (1) stands for the effects of the 
spatially nonuniform Rashba SO interaction, which can 
be written as Hr = H R + V, where H R is the Rashba 
term for A(r, t) = 0, 

H° R = -\o x {V„ A(r)} + l -a y {V x , A(r)} , (4) 

with { , } denoting the anticommutator and A(r) being 
the random Rashba SO interaction. The term V, in turn, 
describes coupling of the electron spin to the external 
field A(r,t) via the Rashba field, 



V 



--A(r) (<r x A y - a y A x ) . 



(5) 



Due to the assumed symmetry with respect to z- 
inversion, the spatially averaged Rashba interaction van- 
ishes, (A(r)} = 0. We assume that the random Rashba 
field can be described by the correlation function related 
to fluctuating density of impurities near the QW ) 24 ' 26 

C XX (r - r') = (A(r) A(r')> = (A 2 ) F (r - r') , (6) 

where the range function F (r — r') depends on the type 
of disorder. We assume the correlator of random Rashba 
interaction in the momentum space in the for m 26 ' 31 
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where R is the spatial scale of the fluctuations. 

In the absence of external field and random Rashba SO 
interaction, the Hamiltonian Hq + Hd describes the spec- 
trum of spin-polarized electrons, e^ a = (fc 2 + fc 2 ) /2m + 



III. NONLINEAR SECOND-ORDER SPIN 
CURRENT 

In the following we consider the z-component of a pure 
spin current flowing along the x axis, that is the only 
component allowed by symmetry of the system under 
consideration. The operators of the electron velocity v x 
and the corresponding spin current tensor component j x 
are 



= i[H 
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where the a-related terms correspond to the anomalous 
contribution to the velocity. The macroscopic spin cur- 
rent density is then given by 



k J 



G 



(11) 



where G^ £ is the Green's function of the system inter- 
acting with the external electromagnetic field. 

Upon substituting (10) into Eq. (11) one can note that 
the second term describes the current caused only by the 
electron density, 



/ — GkF, 

k 



(12) 



conserved under any external perturbation. This conser- 
vation is achieved in calculations by an appropriate shift 
of the chemical potential fi. In equilibrium, however, 
there is no spin current in the system (expected, e.g., 
for QWs with other crystallographic orientations) since 
the integrated contributions from k x /m and a terms in 
Eq. (p~0|) exactly cancel each other. As a result, this type 
of structures does not demonstrate the Rashba paradox 
of the non-zero equilibrium pure spin current^ 2 - This can 
be seen directly by calculating spin current using Eq. (fTTj) 
with the equilibrium Green function in Eq. ([8|l. or by tak- 
ing into account the fact that Hamiltonian in Eq. ^) can 
be transformed by the SU(2) rotation to the form that 
does not have spin dependent terms ^ 
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FIG. 1: Feynman diagrams for the excited spin current in 
Eq. (|lip . The circles correspond to the matrix elements V kk / 
and VW 

However, a nonzero pure spin current, which is the sub- 
ject of interest here, can be generated by an external field 
in the presence of random Rashba coupling, as presented 
schematically by the Feynman graph in Fig. [TJ In this 
graph we introduced the following notations: 

Vide' = A kk / (a x A y (u>) - a y A x (u>)) , (13) 
Vk'k = A k / k (a x Ay(-ui) - a y A x (-uj)) , (14) 



for the transition matrix elements due to the external 
field and spin-orbit coupling. 

With the Feynman graph shown in Fig. [T] we find the 
spin current density omitting the a term in the velocity 
operator since its contribution is conserved if no addi- 
tional electrons are injected to the system. The resulting 
expression for the contribution of a single spin component 
to the total spin current becomes 



jxa — / / -Z— k x I A kk ' | 2 G kecr Gk'e+uia' G kecr , 

(15) 

where A is the vector potential amplitude, making the 
injected current independent on the external field orien- 
tation. The total spin current is given by j x — j x . a= \ — 
ix,o-=-i with j x , a =-i = -jx,a=l- 

Using Eqs. (8) and (15), after rather tedious integra- 
tion of the product of the three Greens functions in the 
complex e plane, one obtains 
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Here we introduced the notation: <5+ = 5 ko - + <5 k <er< and 
<5_ = <5 kcr — 5k> a >. Equation (16) shows that the injec- 
tion of spin current is a coherent effect arising due to 
the change in electron wave function under the resonant 
electromagnetic radiation rather than the injection due 
to the two-photon absorption typical in nonlinear semi- 



conductor optics. 

Since the single electron energy e kcr can be written 
as [(k x + cram) 2 + k 2 ]/2m — ma 2 /2, one can shift the 
chemical potential, /i — > /i + ma 2 /2. Then one can write 
j xa in the form, 
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(17) 



4 




0.2 0.4 0.6 0.: 
0)/|I 



FIG. 2: Total spin current, calculated by using Eq. (|17|l for 
indicated values of /cf-R. The parameters used in numerical 
calculations are given in the main text. 



We emphasize here that the calculated spin current 
is a stationary coherent nonlinear effect proportional to 
the intensity of incident radiation,— in contrast to the 
spin current generated by pulse excitations, where the 
result is proportional to the total fluence in the pulse.— 
The transition processes produce real holes in the initially 
occupied subbands and electrons in those initially empty, 
changing the real occupation of the spin-up and spin- 
down states. The calculated current is also not related to 
the Drude-like linear response at frequency oj, suppressed 
by the factor of the order of (S + /uj) . 

The results of numerical calculation of the spin cur- 
rent (taking part of Eq. (JTTJ) ) are presented in Fig. [5] for 
different values of the parameter kpR. The parameters 
typical for the (110) quantum wells are: 2am/kp = 0.1 
and = 0.1. For the momentum-dependent <5_, we 

assume a typical value, S-/[i = 0.05. Furthermore, we 
used for GaAs: to = 0.067 mo, where mo is the free elec- 
tron mass, Fermi momentum kp = 1.8 x 10 6 cm -1 (corre- 
sponding to electron concentration 5.2 x 10 11 cm -2 ) and 
/j = 18.5 meV. The spin current in Fig. 2 is presented in 
the units of jo, with jo defined as 



. _ 2m 2 ae 2 A2 (X 2 ) 

C 2 7T 3 kp 



(18) 



Taking into account the relation A 2 = (c/uj) 2 E 2 , where 
E is the electric field amplitude, we obtain jo = 
2m 2 a(\ 2 ) (eE/uj) 2 /ir 3 kp, with eE/ui being the ampli- 
tude of the momentum oscillation of a classical electron 
in a periodic electric field. It is interesting to mention 
that the maximum value of E, which still can be consid- 
ered as a perturbation, is determined by eE/uj ~ kp, and, 
therefore, the maximum of jo is of the order of m 2 a (A 2 ) , 
having the physical meaning of the equilibrium spin cur- 
rent induced by the random Rashba spin-orbit coupling. 

To understand better the physical mechanism of 
the nonlinear spin-current generation we consider a 
schematic picture presenting the electron energy bands 
as a function of k x without Rashba random SO interac- 
tion and without external field, see Fig. [31 As we have 
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FIG. 3: Schematic presentation of the light-induced resonant 
formation of spin holes in the energy bands occupied with 
electrons: (a) spin-split energy bands in GaAs (110) quantum 
well; (b) due to the coupling 14. k / (cf. Eq. (|13|l 'l of the spin-up 
and spin-down states, the effective spin (S z ) in each subband 
decreases. 



already mentioned above, the Dressclhaus SO interaction 
leads to spin splitting of the electron states of a free elec- 
tron gas, which results in the energy bands shown in 
Fig. EHa) as a function of k x (for k y = 0). Even though 
the states ktx) of these bands are spin polarized, the spin 
current in equilibrium is exactly zero. This is related to 
the zero current associated with each of the subbands, 
it-,4.> calculated as the flux of electrons in each subband. 
Obviously, vanishing current j a in the subband a means 
that the spin current j% is also zero. Distortions of the 
energy subbands either due to the random Rashba inter- 
action in Eq. (U) or due to the external field in Eq. @ 
do not break the condition j = 0. 

Our calculations, however, showed that nonzero ma- 
trix elements of field-induced spin-flip intersubband tran- 
sitions appear in the presence of random Rashba cou- 
pling. Accordingly, in the noncquilibrium situation the 
electron states in each subband arc a superposition of 
spin up and down states, so that the resulting state |k±) 
(Fig. |3ljb)) has a smaller effective spin. Such mixing of 
kcr) and |kV) states effectively depends on |k — k'| and 
on Iffkcr — ek' a '±uj\, so that the above-mentioned spin mix- 
ing is different at different parts of the dispersion curves 
£k<r- This is shown schematically in Fig.[3jb) for different 
spin subbands. 

Thus, even though in nonequilibrium the current in 
each subband is zero, the associated spin current is 
not zero anymore. For example, in the £kt band more 
up spins flow in +x than — x direction. Correspondingly, 
in the £kj, band more down spins flow in —x direction. 
This results in the net spin-up current in +x direction. 
Obviously, the direction of spin current is related to the 
sign of Dresselhaus coupling constant a. A remarkable 
change of the wave function by a strong electric field 
causes the injected pure spin current of the order of the 
equilibrium spin currents arising as a result of the Rashba 
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paradoxi 25 ! 32 

IV. SUMMARY AND CONCLUSIONS 

To conclude, we have proposed a new effect of the co- 
herent nonlinear generation of a steady pure spin currents 
in GaAs (110) quantum wells by electromagnetic wave. 
The injected spin current is proportional to the intensity 
of the external radiation, strongly depends on the fre- 
quency, and can be injected in the frequency range up 
to the Fermi energy of the two-dimensional electron gas. 
Physical mechanism of the effect is related to the virtual 



spin reorientation of electrons filling the spin subbands 
split by the Dresselhaus interaction in the presence of a 
randomly varying Rashba coupling. The latter may be 
introduced, e.g. by random doping of the quantum well. 
As a result, a 'spin hole' virtually appears in the subband, 
leading to the light-induced spin current. 
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